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Abstract 

The phase diagram of the XX Z chain under the influence of homogeneous 
and staggered magnetic fields is calculated. The model has a rich phase 
structure with three types of phases. A fully ferromagnetic phase, an an- 
tiferromagnetic phase and a massless phase with partial ferromagnetic and 
antiferromagnetic order. This massless phase has its critical fiuctuations 
governed by a conformal field theory with central charge c = 1. When the 
(T^-anisotropy A is zero the model is exactly integrable through a Jordan- 
Wigner fermionization and our results are analytic. For A 7^ our analysis 
are done numerically using lattice sizes up to M = 20. Our results shows 
that for — 1 > A > ^ the staggered magnetic field perturbation is relevant 

and for 1 < A < ^ it is irrelevant. In the region of relevant perturbations 
we show that the massive continuum field theory associated to the model has 
the same mass spectrum as the sine-Gordon model. 



1 Introduction 

The anisotropic Heisenberg model with spin S = 1/2, or XXZ chain, 
is one of the most studied quantum spin system in Statistical Mechanics. 
Since its exact integrability by Yang and Yang|I| this model is considered 
as the classical example of success of the Bethe ansatz. The XXZ chain 
give us the first example of a critical line with critical exponents varying 
continuously with the anisotropy. More recently [^ ^ with the advance of 
the conformal invariance ideas@], the whole operator content of this model 
was obtanied. The critical fluctuations along this critical line is governed by 
a conformal field theory with conformal anomaly c = 1 and moreover their 
currents satisfy an f/(l) Kac- Moody algebra^. 

The effect of an uniform magnetic field in the phase diagram of the XXZ 
chain 1^ ( see Fig. 1 ) is to extend the critical phase over a finite region delim- 
ited by the critical line of Pokrovsky-Talapov (P.T. line in Fig. 1)[0], where 
the chain becomes fully ferromagnetic. Finite-size studies of this extended 
massless phase showed that it is also described by a c = 1 theory with crit- 
ical exponents now depending on the anisotropy and magnetic field. In the 
antiferromagnetic region ( A < — 1 in Fig. 1 ), where the model is massive 
(noncritical), as the magnetic field increases it destroys the massive phase 
entering again in the massless phase. 

Since the uniform magnetic field does not destroy the exact integrability 
of the quantum chain, the eigenspectra in the whole extended massless phase 
is exact solvable. This is not the case when we introduce other perturbations, 
like a thermal perturbation ( which is related to a dimerization in the model 
) or an staggered magnetic field, where exact integrability is lost. Some years 
ago Den Nijs]^, by exploring the connection between the 8- vertex model and 
the Gaussian model, conjectured that a staggered magnetic field could be a 
relevant or irrelevant perturbation depending upon the anisotropy strength. 
Consequently in the region where this last perturbation is relevant we should 
expect an apposite effect as that of the uniform magnetic field, since it will 
bring the system into a massive phase. Motived by this we study in this 
paper the phase diagram and critical properties of the XXZ chain when 
both magnetic fields are present. The XXZ Hamiltonian with these fields is 



given by 

M 



-| ivi 



(1) 
where af , a^ and a^, i = 1,2, ...,M are Pauli matrices attached at the M sites 

of the chain, A is the cr^-anisotropy parameter and h, hg are the uniform 

and staggered fields, respectively. In this paper we assume periodic boundary 

conditions in fffl), and the lattice size M as an even number. 



Using all the machinery coming from finite-size scaling [|TOl and conformal 
invariance[Q, we will calculate the critical lines, exponents as well the masses 
appearing in the underlying massive field theory around relevant perturba- 
tions. 

The paper is organized as follows. In Section 2 we review the conformal 
invariance relations relevant for our purposes and the operator content in 
the absence of magnetic fields. In Section 3 we consider the cases where 
hg 7^ 0. Initially we consider the case where A = or the XY chain. This 
case is special since the eigenspectra can be calculated analytically through 
a Jordan- Wigner fermionization of the Hamiltonian. The general situation 
where A 7^ is studied numerically. In Section 4 we conclude with a general 
discussion of our results. 

2 Conformal invariance relations and the op- 
erator content when /i^ = 

Like most statistical mechanics systems in 1 -|- 1 dimensions we assume 
that the Hamiltonian (|l|) in its massless regime is conformally invariant. 



Under this assumption for each operator [g, ^ Oa with dimension Xa and 
spin Sa in the operator algebra of the infinite system, there exists an infinite 
tower of states in the quantum Hamiltonian, in a periodic chain of M sites, 
whose energy and momentum as M — ;> oo are given by 

077-7; 

El^,{M) = Eo{M) + -—(x^+j+f) + o{M-') (2) 



and 



2tt 
^^/(M) = -(..+ J-/); (3) 



where j,j' = 0, 1, ..., Eq{M) is the ground-state energy and v is the velocity 
of sound that can be determined by the dispersion relation of the spectra or 
by the difference among energy levels belonging to the same conformal tower. 
The conformal anomaly c can be obtained from the finite-size corrections 
of the ground-state energy. For periodic chains, the ground-state energy 
behaves asymptotically as[T^ 



M °° 6M2 



+ o{M'') (4) 



where Coo is the ground-state energy per site in the bulk limit. 

The above relations has been applied successfully to a large number of 
statistical mechanics systems. In the case of the XX Z Hamiltonian with no 
external fields, H{A,h = 0,hs = 0) and —1 < A < 1 these relations give 
us the whole operator content of the underlying field theory]^, ^ |^. Along 
the critical line — 1 < A < 1 the central charge is c = 1 and the anomalous 
dimensions appearing in the model are given by integer numbers or by 

2 "^^ / \ 

where Xp = ^~'^%7r ^^^ '^ ^^^ "^ ^^^ integers. The dimension Xo,i = Xp 

is the same as that of the polarization operator in the six vertex model. 
Moreover we can show0, by combining the integer dimensions with Eq. (|^), 
the model is governed by a conformal field theory satisfying a larger algebra 
than the Virasoro conformal one, namely, a U{1) Kac-Moody algebra. This 
[/(I) symmetry comes from the commutation of the Hamiltonian (|I]) with 



the z component of the total spin[^ 



^ M 

S. = ^Y.^t = n. (6) 

i 

For a given sector of the Hilbert space where the total spin Sz = n it will 
corresponds a set of primary operators On,m with dimensions Xn,m (^^ = 
0,±1,±2,...) and the number of descendants operators (with dimensions 
Xn,m + j + f 5 jif £ Z) will be given by the product of two U{1) Kac- 
Moody characters^. The operators On,m can be interpreted, in a Gaussian 
language ||14||, as the operators with vorticity n and spin wave excitation num- 



ber m. The critical exponents are obtained from the dimensions (W. For 



example the exponents rjx and rjz governing the correlations 

< (T"(0)a"(r) > ~ r-^- (7) 

< a'{0)a'{r) > ~ r"^^ (8) 

are related to the lowest eigenenergies in the sector with S^ = 1 and S^ = 0, 
respectively. The exponents are given by rj^ = 2a;o,i = ""~^°'^ — ^——^ and r]^ = 
?7~^ for > A > — 1 and rj^ = 2 for 1 > A > 0, since in this last region of 
anisotropics the lowest dimension above zero is one. 

All the above results were obtained for hs = h = and by explor- 
ing the exact integrability of the model through the Bethe ansatz, which 
numerically [0] or analytically]^] produced very precise results. Such inte- 
grability is not destroyed by the introduction of the uniform magnetic field 
and precise results can also be obtained in this case. In the case where 
hs = 0, h ^ 0, as we see in Fig. 1, the massless line extends forming a 
massless phase. This massless phase is ordered ferromagnetically and for 
a given value of the anisotropy A the magnetization changes from zero to 
its maximum value where the system is fully ordered, the points where the 
magnetization reaches its maximum value form the Pokrovsky-Talapov tran- 
sition line[0], which separate the massive fully ferromagnetic phase from the 
massless partially ordered ferromagnetic phase. 

The physical mechanism producing the extended massless phase in Fig. 
1 is due to the following. From Eqs. (|I]) and (^) the magnetic field only 
decreases the eigenenergies in the sectors where n > 0. In the critical line 
— 1 < A < 1 and h = the ground-state belongs to the sector n = 0. 
Since we have no gap a small but finite positive magnetic field will lead 
the ground-state into a sector with n > 0, producing a partially ordered 
ferromagnet phase. However the fluctuations around this ordered phase are 
of the same nature as those appearing in the absence of the magnetic field 
{h = 0). This argument imply that we should expect in the whole massless 
phase a c = 1 Gaussian field theory with the anomalous dimensions like 
those given in Eq. (|^). In fact this fact was observed analytically §] for 



— 1 < A < 1 and numerically! 15] in the whole extended phase. In Refs. [g], 
although a closed analytic form for the dimensions in the extended massless 
phase was not derived, the dimensions are shown to be of Gaussian type like 
in (1^), but with Xp now depending on A and h. The case where A = is 
special since only in this case the dimensions do not depend on the magnetic 
field h. 



The finite-size analysis of tlie Hamiltonian (|I|) witli a nonzero uniform 
magnetic field deserves some comments since for a fixed value of h tlie ground- 
state clianges sector as tlie lattice size increase. In order to keep the bulk limit 
physics (M -^ oo) we should keep the magnetization per particle fixed in the 
finite chains. This imply that for a given magnetic field only a sequence of 
lattice sizes producing the fixed magnetization, related to /i, should be used. 

The transition line of Pokrovsky-Talapov type in Fig. 1 is obtained by 
calculating the lowest value of the magnetic field where the ground-state is 
the ferromagnetic state ( sector n = M/2 ). This magnetic field is /i = 1 — A 
for arbitrary values of the lattice size M. 

In the following sections we will calculate the effects where the staggered 
and the uniform magnetic field are included. 

3 Results for hs ^ 

We now analyze the situation where both magnetic fields h, hg are 
nonzero. The staggered magnetic field for h = produces an antiferro- 
magnetic ordered ground-state. The ground-state belongs, as when hg = 0, 
to the sector with n = 0. In this section we are going to study analytically 
for A = 0, and numerically for arbitrary A the nature of the resulting phase 
when hg is introduced. Irrespective of the hg value we do expect that for 
sufficiently high values oi h {h » hg ) the same mechanism discussed in last 
section which changes the ground-state sector will also take place. Conse- 
quently we should have a phase transition surface of Pokrovsky-Talapov type 
where the ground-state becomes the fully ordered ferromagnetic state. The 
magnetic field h{A, hg ) at this surface, as in the previous section, is calcu- 
lated by imposing the ground-state as the fully ordered ferromagnet state. 
This critical surface 

h = hr = ,/}^Ti-A. (9) 

is obtained by equating the lowest eigenenergy Ei = —^{M — 4) — h(M — 
2) — '2, J hi -|- 1, in the sector with one spin down with the ferromagnetic 
energy (zero spin down) Eq = — (^ + h)M. For h > h^^ the ground-state 
of H{A, h, hg ) is the ferromagnetic state. In order to study the cases where 
h < h^'^ we will consider first the special case where A = 0. 



3.1 A = 

In this case the Hainiltonian H{A, h, hg ) can be diagonahzed through 
a Jordan- Wigner transformation |T^. This is done by introducing fermionic 
operators 

cj = II e^""' "' a- , c] = a+ [] e"*""' "' (10) 

1=1 1=1 

where af = ^{a'j ± iaj) for j = 1,2, ...,M. It is easy to check that these 
operators satisfy the fermionic algebra 

{c], Cy} = 5jj' {c], c],} = {Cj, Cy} = . (11) 

In terms of these operators the conserved U{1) charge is given by 

1 ^ ^^ M ^^ M 

^. = ^E^r = E--^-r-y = E4c.-f (12) 

i 1=1 i=l 

which imply that the Hilbert space associated to if (0, h, hs ) can be separated 
into block disjoint sectors labelled according to the fermionic number 

n^=y: cU = — ^ ■ (13) 

For the sake of simplicity let us restrict to the case where M is multiple of 
4, so that Np and n have the same parity. For a given number of fermions 
Nf the Hamiltonian is given by 

Af-l 



H{0,h,h) = Y. [ -(4q+i + c1+iq) + (-1)^/1, (2cJq-1) 

-2n/i + e^'^^^ (4ci + cIcm) • (14) 

We see that the sectors with even and odd number of fermions are obtained 
by antiperiodic and periodic boundary conditions respectively. 

Since the Hamiltonian (0) is quadratic in fermions operators we can (see 
Appendix A of ref. [0), through a linear combination of Cj,Cj introduce 
new fermionic momentum operators f]k,Vk ^^^^ brings the Hamiltonian into 
a diagonal form 



H{0, h,hs)=J2 ^dkJhl + cos2 k r]lr]k - 2nh (15) 



Ok = +1 for 



71 . n 

<k< - 

2 - 2 



and 6k = —1 



otherwize 



(16) 



where the sum now runs over the momenta —it < k < it and 6^ are signals [0]. 
Due to the boundary condition term in (11^) the possible values of momentum 
k will depend on the lattice size parity and are given by 



(2' + ^i- 



for Np even, and 



/ 



M 

Y 



M 

y 



- — ,-— + 1,...,— -2,-^ 



M 

Y 



M 

Y 



(17) 



, TT , MM M M 

21—, 1 = , + 1, ..., 1,— 

M ' 2 ' 2 ' ' 2 '2 



(18) 



for Np odd. The eigenenergies in a given sector with Np fermions are ob- 
tained from the jy (mIn v. combinations of fermion energies. The momentum 
of a given eigenstate is obtained by adding the fermions momenta P = X]fc k. 
For example the lowest eigenenergy in the sector with U{1) charge n will be 



Er, 



for n even and 



E„ 



(2' + i'¥ 




h^ + cos2 ( 21 — 



2nh 



-2nh 



(19) 



hi 



cos' 



'- 4+2+2 



21- 



TT 
'M 



-2nh 



-2nh, 



(20) 



for n odd. These energies corresponds to zero-momentum states since they 
come from a symmetric distribution of fermion momenta. 
From 



^^ we know that at h = Jhl + 1 the system undergoes a phase 
transition to the ordered ferromagnetic state. In order to detect other phase 
transitions let us calculate the mass gap of if (0, hg , h) for small values of 



h. In such cases the ground-state remains in the sector n = and the first 
excited state in the sector n = 1. From Eqs. (0) and (^) the mass gap is 
given by 

G = hm {El - Eo) = 2(/i, - h) . (21) 

Therefore, as long as hg > h the model is massive and ordered antiferro- 
magnetically. At h = hg the model becomes massless and by the same 
mechanism we discussed in the last section we expect that as h increases the 
model stay massless until the Pokrovsky-Talapov line is reached. This inter- 
mediate phase has a ferromagnetic order since now the ground-state change 
sectors as h increases. In Fig. 2 we show the phase diagram at A = 0. 

Intermediate phase at A = 

In order to analyze this phase using finite-size scaling we should keep the 
magnetization per particle fi fixed. This is done by choosing a sequence of 
lattice sizes Mi,M2, ... and sectors ni,n2, ..., such that fi = ff- , i = 1,2,... 
is fixed. For a given lattice size Mi and magnetization /i the magnetic field 
that brings the ground-state in the sector where n, = /xMj should be in the 
range hminiMi, fx) <h< h^ax{Mi,^i), where hmin{Mi, fx) = (E„^ - E„,^_i)/2 
and hmax{Mi, fi) = {En.+i — En-)/2. From Eqs. (|TUp and ( PD| ) we obtain in 
the bulk limit 



hmin{00, /i) = hraax{00, fi) = h^ = ^ k^ + slu a , (22) 

where f^ = j^ and a = fin. Consequently if the applied field has the value 
h = h^ all the finite lattice sequence Mj (z = 1,2,...) will have its ground-state 
in the sector with U{1) charge rti = fiMi {i = 1, 2, ...). It is simple to verify 
form Eqs. (|l9|) , (pOD that for arbitrary values of < ^u < 1 the gap vanishes as 
M — i> oo, and consequently this partially ordered phase is critical. Assuming 
the critical fluctuations in this phase governed by a conformal field theory, 
the relations (@),(0) and (^ of Sec. 2 enable us to calculate the anomalous 
dimensions and central charge of this underline field theory. 

For fixed value of /i the finite-size correction of the ground-state energy 
E^^j in the asymptotic regime M —>■ oo can be obtained by using the Euler- 
Maclaurin formula 1 18 1 in Eq. (|T9|), which give us 



^?M TT sin(2a) f ..-2\ /oq\ 



M QM^ h^ 



where 

Coo = / Jhl + cos2 1 dt - 2hafi (24) 

is the energy per particle in the bulk limit M ^ oo. The low lying excited 
states with nonzero momentum are calculated by changing the fermions mo- 
menta entering in the ground state. From these energies we can derive the 
dispersion relation and the sound velocity 

C=5lf^ (25) 

which changes continuously along the massless phase. By comparing Eqs. 
( p3D and (|25|) with Eq. (|^) we see that the conformal anomaly is c = 1 for 
all values of fi. 

From Eq. (^ the anomalous dimensions of the conformal operators are 
calculated from the finite-size corrections of the mass gap amplitudes. Such 
corrections, for a certain fermionic configuration, are calculated straightfor- 
wardly by using the Euler-Maclaurin formula. For an arbitrary but rational 



magnetization /z = 2 (^or magnetic field h = Jhl + sin /z) the lattice se- 
quence with size Mj = iq {i = 1, 2, ...) will have its ground-state in the f/(l) 
sector Hi = fiMi = ip (i = 1,2, ...). The lowest eigenstate in neighboring 
sectors where nl = Ui + n {n = ±1, ±2, ...) has zero momenta and will give 
the mass gap 

The excited states in the sectors n^ {i = 1,2,...), obtained by the addition 
(subtraction) of a momenta j^ if^] ^o the fermions with positive (nega- 
tive) momenta will give us an eigenstate with momentum 

(27) 







pi^f 




and mass gap 








E^f 


rpGS 




( n' 



(m-1 



(28) 



The relations @ and (||) enable us to identify the eigenenergies E^''' as the 



1^2 



conformal tower of the primary operator with dimension Xnfl = in 

9 



Other excited states in the sectors where n^ = rii -\- n {n = ±1, ±2, ...) 
are derived from the fermion configuration producing i?;^';' . These eigenen- 
ergies are obtained by replacing fermions with positive (negative) momenta 
by fermions with negative (positive) momenta. The momenta of such states 
are 

P^;^ = 2nix+'^inm + j-f) (29) 

and the mass gaps are 

Bii - ESf = |C f n% +^+J+f)+o (M-') (30) 



P 



where Xp and n,m,j,j' = 0, ±1, ±2, ... . From relations (|]) and (||) we iden- 
tify Eq. (|30|) as the conformal tower of the primary operator with dimension 
Xn,m = (x + "^^)- Following the discussion of Sec. 2 ( see Eq. (||) ) we can 
interpret Xn,m, as the dimension of a Gaussian operator On,m with vorticity 
n and spin-wave number m. Furthermore, from the degeneracies of the low 
lying energies in the conformal towers \ E^^m \ we verify that their are given 
by the product of the characters of f/(l) Kac-Moody algebra as discussed 
in Sec. 2 for the case where hg =0. Therefore the whole massless regime 
has the same critical nature as in the case where hg =0, being governed by 
a c = 1 conformal field theory with operators satisfying a U{1) Kac-Moody 
algebra. 

Comparing Eqs. (30), (25) and (|^) we see that while the sound velocity 
depend on h and hg the dimensions Xn,m are constants. For example the 
correlations (|^ and (|]) are given by 

< a^(0)a^(r) > ~ r'l (31) 

< (T'{0)a'{r) > -fi^ ~ r-2 (32) 

where the magnetization per site /i is given by 

/i = - sin-i (v/^' - hA . (33) 

As we shall see this independence of exponents is a particular feature of 
the case where A = 0. It is interesting to observe that contrary to stan- 
dard conformal towers appearing in critical systems, Eq. (129|) tell us that 

10 



the eigenstates in the conformal tower has a macroscopic momenta as long 
/i 7^ 0. This is due to the fact that the fluctuations ruling the critical be- 
havior are in top of a partially ordered ferromagnet state. Another delicate 
point appearing here[|19| is that a continuum conformal theory can only be 



construct for h = J hi + sin^ /i, with /i having a rational value, since the lat- 
tice sizes and sectors appearing in the finite-size sequences are integer values 
fJ' = jj- (i = 1,2, ...). For h connected with an irrational value of fj, we 
can only obtain a conformal theory by approximating /i by a close rational 
number. 

Massive phase at A = 

The antiferromagnetic phase for h < hg is massive with a mass gap 
given by (|2l|) . A continuum field theory describing the physics in this mas- 
sive phase can be obtained in the neighborhood of the perturbing parameter 
6 = hs — h >0. Such field theory will be massive and the masses can be 
estimated from the finite-size behavior of the eigenspectra. The mass spec- 
trum can be calculated by applying the scheme followed by Sagdeev and 
ZamolodchikovpOH in the study of the Ising model in an external magnetic 



field. To do such calculations we should initially find the finite-size cor- 
rections of the zero-momenta eigenenergies Ek{6,M), k = 1,2,..., at the 
conformal invariant point 6 = 0. In the case of the Hamiltonian (|I]) the 
results of Ref. P] tell us that such corrections, not only at A = but for 
arbitrary values of A (—1 < A < 1), are governed mainly by the irrelevant 
operator with dimension x = xo,2 = 1/a^p and the descendent of the identity 
operator with dimension 4. From Ref. [0 we have 

Ek{6 = 0,M) = e^M + —-{xk-—)+a^{ — ) + (34) 



M V 12 / \MJ 

1 \ x—l / 1 \ 2S— 3 / 1 \ 3S— 5 



where Xk is one of the dimensions (BI) associated to E^ and ai, 02, ... are M- 



independent factors. According to the scheme of Ref. ||20[ if the perturbed 
operator which produces the massive behavior has dimension y we should 
calculate the eigenspectra in the asymptotic regime 5 — > , M — > 00, with 

X = 6^ M (35) 

11 



kept fixed. In this regime (53) is replaced by 



Eki5,M) = e^M + 5—F,iX) + S 



2^(--i) 



2-y 



GkiX) + 



(36) 



52^(2--3)y^(^) + S2h(^^-^')HkiX) + ... . 



The masses in the massive continuum field theory are obtained from the 
large-X behavior of the function |Q Ffc(X), i.e., 

ruk ~ Fk{X) - Fo{X) (37) 



The relations in Eqs. (p^)-(|37|) can be used for arbitrary A. In fact once 
we know x = 1/xp, the dimension associated to the off-critical perturbation 
y as well the generated masses can be calculated from Eqs. (p5|)-(|37|). 

In the case A = the perturbing parameter is 6 = hg —h >0. By keeping 

X = [hg — h)^^ M fixed, with y unknown, we obtain from Eqs. (|l^) and 



(3^ 



Ei-Eo = 2{K -h)- 



2 E 

/=1 



Ki + sin^ 



27r ^- 



Ki + sin^ 



TT 



(21 - \)h 



2-!/ 



which gives for hg 







E^-E^^2{hg -h)-%hr- 



(39) 



where fi is a constant. Consequently by comparing Eq. ( PU] ) with Eqs. 

( |37D we obtain the dimension y = 1 for the perturbing operator and the mass 

mi = 2 associated to the first gap. 

The gap associated to the lowest eigenenergy in the sectors with n even 
is obtained from Eq. (|19D 



n/2 \ 



TT 

X 



(2/ - \)h. 



1 1 

2-!/ 
S 



(40) 



Since n is finite, expanding for small values of hg we get, by comparing with 
Eq. (l36|), again the dimension y = 1 for the perturbing operator and 



n/2 

F„(X)-Fo(X)=4^ ^/l 



TT 

X 



(2/ 



(41) 



12 



which gives, for X — > cxd, the masses rUn = n 2 = n mi. The same result 
is also obtained in the case where n is odd. Similar calculations for zero- 
momentum (modulo tt) excited states show us that these states are associated 
with multiples of the mass mi. However our calculations shows the existence 
of two equal masses in the system : mi = m2 = 2. This can be seen clearly by 
a two-fold degeneracy of the first excited state with zero momentum (modulo 
tt) in the ground state sector n = 0. The levels are related to the threshold 
energy 2m i. 

The dimension y = 1 for the perturbing operator and the Gaussian di- 
mensions given in Eq. (^, with Xp = 1/4, indicate that this operator is the 
Oo,i operator (see Sec. 2) with dimension y = xo,i = 1. This imply that in 
the region where a;o,i > 2, (A < ^) this perturbation will be irrelevant and 
the phase stays massless after the introduction of a small hg . 

3.2 A 7^0 

In this case the analytical diagonalization done for A = does not work 
since the Hamiltonian through the Jordan- Wigner transformation (|I0|) will 
have a four-body fermionic interaction. Consequently our analysis will be 
done numerically. The critical surfaces of H{A,h,hs ), according to finite- 
size scaling theory (FSS) [|1^ can be obtained from the extrapolation (M — »• 
oo) of the sequences (A^*^), h^^^\ h'^^^) (M = 2,4, ...), obtained by solving 
the equation 

MGm(A, h,K) = {M- 2)Ga/-2(A, /z, hs ) (42) 

where Ga/(A, h, hg ) is the gap of the Hamiltonian ([l|) with M sites. 

In order to simplify our analysis let us consider initially the plane h = 0. 
The critical surface is obtained by solving (^21) with A or hg fixed. In Table 1 
and 2 we show some of the finite-size sequences for lattice sizes up to M = 18. 
In Table 1 hg is fixed while in Table 2 A is fixed. In Fig. 3 we show the 
curves in the space (A, /i), that solve Eq. (|^) for several lattice sizes. We 
clearly see in the figure that as (M -^ oo) the curves for A < A* ^ 0.7 tend 
toward hg =0, while for A > A* the curves tend toward nonzero values 
of hg . The blow up region A f» A* is also included in Fig. 3, where the 
points we used to draw the continuum curves are also shown. These results 
clearly indicate that for A < A* ^ 0.7 the perturbation introduced by hg is 
relevant producing a massive antiferromagnetic phase, while for 1 > A > A* 

13 



the model stays massless until it reaches a critical staggered field /i* = h*{A) 
which depends on the value of A. This is in complete agreement with the 
predictions we did considering the A = results, where we associated the 
operator with dimension xo,i = l/xp = 2(7r-cos-i(-A)) ^i^h the staggered 
field perturbation. Therefore the perturbation is irrelevant for A > A* = 
^ ~ 0.707 and the phase which appears for A > A* in Fig. 3 should be a 
massless antiferromagnetic ordered phase, since hg 7^ in this phase. This 
is complete agreement with predictions did in an earlier work by den Nijs]^]. 
The massless phase (denoted by ML in Fig. 3) is separated from the fully 
ferromagnetic phase (denoted by FE in Fig. 3) by the Pokrovsky-Talapov 
line (PT) hs = VA^ - 1, obtained from Eq. (1). 

We confirmed numerically the massless nature of the ML phase. Using 
relations (@),(0) and (|) we verified that the whole phase is governed by a 
c = 1 Gaussian conformal field theory with dimensions as in (|]), with Xp 
varying continuously in the whole phase. As a consequence of (^), ratio of 
dimensions like Xnfl/xifi = n^ should be independent of Xp. In order to 
illustrate this let us consider the ratio G2IG1 between the gaps in the sectors 
n = 2 and n = 1. Since the gaps are related to the dimensions X2,o and Xi^, 
as (M -^ cxd) this ratio should tend to the value X2,o/a^i,o = 4. In Fig. 4 
we show these ratios as a function of A for M = 20 and some values of hs . 
It is clear from this figure that the ratio changes to the expected value of 4 
at values of A which depend on the hg value, which is consistent with the 
results shown in Fig. 3. 

Our results also show that the PT line of Fig. 3 has a different behavior 
than the rest of the Pokrovsky-Talapov surface (^, since only along this line 
we have a first order phase transition. When h = hj crossing the PT 
line the ground-state changes from the ferromagnetic sector n = M/2 to the 
sector n = discontinuously. As a consequence the magnetization per lattice 
point changes discontinuously from 1 to zero, as we move from the FE phase 
to the ML phase (see Fig. 3). Similarly the staggered magnetization 



(s:) = o E i-y^! (43) 




changes discontinuously from zero to a finite value, which depend on A. Also 
along the surface PT the ground-state wave function which is nondegenerated 
for h ^ becomes M-degenerated when h = 0. All the lowest states in the 
sectors with n = — ^,— ^ + 1,...,^ — 1,^ become degenerated. We can also 
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show that along this PT curve these degenerated ground-state wave function, 
for a given sector n, is given exactly by 

^„(A,/i, ,/i = 0)=^'g3S"i(-)'^> \si,S2,...,sm) (44) 

where 

A = ^/^^H^=^, (45) 

Si = ±1 {i = 1,2,..., M) are the eigenvalues in the az basis and the prime in 
the sum indicates that the configurations {s} are constrained to J^tii Si = n. 

Massive phase at A 7^ 

Similarly as we did for A = let us now investigate the continuum massive 
field theory describing the fluctuations in the massive phase (AF) of Fig. 3. 
The masses of this continuum theory are obtained from the eigenspectrum 
of the Hamiltonian in the scaling regime given by Eq. (|35|) with 6 = hg 
and from our present analysis the dimension y of the perturbing operator is 
y = a;o,i = 2(7r-cos-i(-A)) - Using Eqs. (|3|), (H) and {^ we can calculate 



the mass ratios of the continuum theory. These are calculated from the 
asymptotic regime M —>■ 00 and X ^ 00 of the finite-size sequences 



,{«) 



Fi"^(X,M)-Fo^°^(X,M) 

Fi^^ (X, M) - Fo^°^ (X, M) ml 



Rrix,M) = ^'j^' Jo) J„; ^ 7:J (46) 



The functions F^ (X, M) are obtained by using in (|36| ) the finite-size se- 
quence of the k zero-momentum state {k = 0, 1,2, ...) in the U{1) sector n. 
The associated mass is m}^ while the lightest mass, obtained from the lowest 
eigenenergy in the sector with n = ±1, is mo • Our results indicate that for 
A > we have only two equal masses M= itiq = rriQ and a continuum 
starting at 2M. For A < other masses appears depending on the value of A. 
In Table 3 we present some of our estimators for the mass ratios. We show 
for M = 20 the ratios ( PB| ) for some values of X and A. Our results show the 
appearance for A = —1/2 of two masses mi = mo and m2 = mi with a 
ratio m2/mi = 1.61 ± 0.03. For A = — ^ (A = — ^) a third mass appears 

with the ratios m^/rrii = ttiq /itiq = 2.0 ± 0.2 (= 1.9 ± 0.1) ,m2/mi = 
1.41 ±0.01 (= 1.2 ±0.1). 
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It is well know that the six-vertex model (or the XX Z chain) is trans- 
formed into the eight-vertex model (or the XYZ chain) by a thermal 



perturbation|21]. The spectrum of the transfer matriz of the eight-vertex 



model was calculated |2^ explicitly and is the same as that of the sine-Gordon 
model [^. From these results the thermal perturbation in the XX Z chain 
will produce the following sine-Gordon masses. For A > there will have 
a mass M twice degenerated plus a continuum starting at 2M. These are 
the soliton-antisoliton pair. For A < 0, beyond the soliton-antisoliton pair 
bounded states will appear with masses depending on A, 

mi+i = 2mi sinf --j , i = 1, 2, ..., [a] (47) 

where 

37r-4cos-i(-A) ,_,. 

a = ^ ^ (48) 

7r 

and [a] is the integer part of a. 

Our results of Table 3 shows that also in the case where the perturbing 
field is the staggered magnetic field the mass spectrum is the same as the 
sine-Gordon model given by Eqs. (E^) and 



4 Conclusions 

The results we obtained in last section enable us to obtain a complete 
schematic phase diagram of the Hamiltonian (|I]). For simplicity we only draw 
the phase diagram for A > in Fig. 5. The surface OGFHI separates the 
antiferromagnetic phase AF from the massless phase ML , while the surface 
AEDCB separates this last phase from the frozen ferromagnetic phase FE. 
The order parameter characterizing such phases is the magnetization 




Mk = { -E^r ) (49) 

and the staggered magnetization 

M,^ = ( ^Y.i-y^t ) ■ (50) 
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In phase FE M,, 7^ and Mh, = 0, in phase AF M/^ = and 
'Mhs ¥" 0, while for /i 7^ in phase ML we have M/x 7^ and 'Mhs 7^ 0. At /i = 
phases AF and ML has only antiferromagnetic order {.Mh = 0, Mh^ 7^ Oj, 
but these phases differs since contrary to phase AF, (or FE) where the corre- 
lations have an exponential decay in the massless phase ML the correlations 
has a power-law decay. 

All the phase transition surfaces of Fig. 5 are of second order except 
the curve AE obtained when /i = 0. The surface AEDCB is the Pokrovsky- 
Talapov surface given by Eq. (^. In this whole surface the ground state 
is the ordered Heisenberg state with all the spin up, except at the line AE 
where the ground-state is M degenerated with a wave function given exactly 
by Eq. (HD. 

In the entire massless phase our numerical and analytical results show 
that the critical exponents change continuously with /i, hs and A, except at 
A = where the exponents are constant. The critical fluctuations in this 
phase is governed by a conformal field theory with central charge c = 1 and 
operators satisfying a U{1) Kac-Moody algebra. 

Our numerical and analytical analysis clearly indicate that the staggered 
magnetic field perturbation is associated to the operator with dimension 
2^0,1 = 4^ = 2(7r-cos-i(-A)) ' '^^^^ confirms an earlier prediction of den Nijs^] 
which states that for A > ^ this perturbation is irrelevant. 

Coming from the massless phase ML to the antiferromagetic phase AF, 
of Fig. 3 the fluctuations will be ruled by a massive continuum field theory. 
Our analysis of last section indicate that such massive field theory is the 
sine-Gordon field theory with a mass spectrum that depends on the value of 
A like in Eq. (^Tf). This is the same mass spectrum obtained by perturbing 
the XX Z by a thermal field. 
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Figure Captions 



Figure 1 - Phase diagram of the Haniiltonian H{A, h, hg ) given in Eq. 
(|l]) for hg = 0. The different phases and the Pokrovsky-Talapov hne are 
indicated. 



Figure 2 - Phase diagram of the Hamiltonian H{A, h, hs ) given in Eq. 
(|I]) for A = 0. The massless phase is separated from the ferromagnetic phase 
by the Pokrovsky-Talapov hne h = J hi + 1 and from the antiferromagnetic 
phase by the hne h = hg . 



Figure 3 - Finite-size estimators of the phase diagram of the Hamiltonian 
if (A, /i, hs ) given in Eq. (|l]) for /i = 0. The line are obtained by solving 
Eq. (gg) for lattice size pairs (M - 2, M) (M = 10, 12, ..., 18). The antiferro- 
magnetic, massless and ferromagnetic phases are indicated by AF, ML and 
FE respectively. The Pokrovsky-Talapov line PT is also indicated. The blow 
up region 0.57 < A < 0.82 is also shown with the points used to draw the 
continuum curves. 



Figure 4 - Ratios G2/G1 between the gaps G2 in sector n = 2 and Gi in 
the sector n = 1 of the Hamiltonian H{A, 0, hs ) given in Eq. (|1[). The ratios 
are shown as a function of A for lattice size M = 20 and some values of hs . 
In the region where the phase is massless we expect the Gaussian relation 

G2/G1 — ^ 3^2,0/3^1,0 = 4 as M — ^ 00. 



Figure 5 - Schematic phase diagram of H{A, h, hs ) given in Eq. (|1]). The 
phases indicated by FE, ML and AF are the fully ferromagnetic, massless 
and antiferromagnetic phases, respectively. The Pokrovsky-Talapov surface 
ABODE separates the ferromagnetic and massless phases and the surface 
FGOIH separates this last phase from the antiferromagnetic one. All the 
phase transition surfaces are of second order excepting the line AE where 
the transition is first order. 
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Table Captions 



Table 1 - Sequences of estimators for the anisotropy A in the plane h = 0, 
obtained by solving Eq. p^) with hs kept fixed. 



Table 2 - Sequences of estimators for the anisotropy A in the plane h = 0, 
obtained by solving Eq. (|4^) with A kept fixed. 



Table 3 - The mass-ratio estimators R}^ {X, M) defined in Eq. (^Gf ) for 
some values of the anisotropy A in the plane h = 0. The conjectured values 
of last line is given by Eq. (^^ . 
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Table 1 



M-2,M 


hs = 0.01 


hs = 0.03 


hs = 0.05 


hs = 0.1 


hs = 0.2 


hs = 0.5 


hs = 1.0 


8,10 


-0.4540731 


-0.0159442 


0.3100008 


0.5969325 


0.7614057 


0.9864091 


1.3628443 


10,12 


-0.3507805 


0.2311949 


0.4642322 


0.6591962 


0.7855170 


0.9930197 


1.3634265 


12,14 


-0.2052818 


0.3724661 


0.5441739 


0.6931414 


0.7996157 


0.9972828 


1.3660457 


14,16 


-0.0457780 


0.4570385 


0.5922408 


0.7145697 


0.8088957 


1.0002445 


1.3650563 


16,18 


0.0891170 


0.5122369 


0.6242583 


0.7293916 


0.8154915 


1.0024008 


1.3657592 



Table 2 



M-2,M 


A = 0.6 


A = 0.65 


A = 0.675 


A = 0.7 


A = 0.725 


A = 0.75 


A = 0.8 


8,10 


0.1010376 


0.1213905 


0.1344105 


0.1499071 


0.1682486 


0.1890890 


0.2420444 


10,12 


0.0776255 


0.0958401 


0.1078859 


0.1225705 


0.1404855 


0.1621757 


0.2176779 


12,14 


0.0624300 


0.0789604 


0.0901779 


0.1041465 


0.1215976 


0.1432584 


0.2005065 


14,16 


0.0518302 


0.0669150 


0.0773753 


0.0906412 


0.1075649 


0.1290517 


0.1875885 


16,18 


0.0440532 


0.0578975 


0.0676763 


0.0802761 


0.0966522 


0.1178787 


0.1773956 



22 



Table 3 





A = 4 


A = 


V2 
2 


A = 


2 


X 




7721 


"la 
mi 


mc; 
mi 


m.s 
mi 


6 


1.64115 


1.41970 


2.20667 


1.25317 


2.15167 


10 


1.61967 


1.40502 


2.06304 


1.18785 


2.01352 


14 


1.62223 


1.40555 


2.03037 


1.15271 


1.99746 


Eq.® 


1.61803 


1.41421 


2.00000 


1.24698 


1.94986 
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